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Особливості трансформацій гравітаційного поля в задачі Алексідзе 
Результатом зйомок гравіметричного поля є значення 2(5) модуля градієнта потенціалу сили тяжін- 


ня (МГПСТУ) на деякій гладкій замкнутій поверхні 5, яка цілком охоплює всі тяжіючі маси. Напрям 2(5) 


невідомий, тому задачу трансформацій поля неможливо звести до лінійних граничних задач. Наближений 
однозначний розв'язок вказаної задачі можна отримати і на незамкнутій проте достатньо великій (поряд- 
ку 19х1?) поверхні 5. 

Спосіб трансформації. Один із шляхів трансформації і перерахунку значень сили тяжіння пов'язу- 
ють з розв'язанням нелінійної зовнішньої граничної задачі Алексідзе для потенціалу сили тяжіння || Ї: 
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Цю задачу лінеаризують шляхом певних припущень щодо напрямку градієнта потенціалу на поверхні 5. 
Такі припущення вносять певні похибки, бо визначення напряму сили тяжіння - громіздка задача, чи не 
складніша за визначення модуля сили тяжіння. 

Проблема розв'язності і єдиності цієї нелінійної граничної задачі (1) принципова для гравіметрії: 
відповідає, чи можливе точне однозначне визначення сили тяжіння в зовнішньому просторі на основі то- 
чного визначення абсолютного значення сили тяжіння на замкнутій поверхні 5, що цілком охоплює усі 
маси і яка додаткова інформація потрібна для цього. За одночасного визначення значення і напряму сили 
тяжіння на поверхні 5 обчислення сили тяжіння в зовнішньому просторі є розв'язком задачі Діріхле. 
Сумісні визначення модуля і напряму сили тяжіння дорогі 1, в масовій гравірозвідці не проводяться. 

Розв'язують граничну задачу (1) методом послідовних наближень. У плоскому випадку задача (1) ви- 
значення потенціалу, як гармонічної функції, модуль градієнта якої заданий на замкнутому контурі, не 
має єдиного розв'язку. Якщо дійсна частина и(х,у) голоморфної на сфері С, яка повністю охоплює 5, 
функції Д(с)зи (ху) З і» (х,у) є розв'язком граничної задачі, уявна частина теж є розв'язком тієї ж гра- 
ничної задачі. 

Для розв'язання граничної задачі (1) у плоскому випадку задіяно такий ітераційний процес 
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дей" - К-те наближення розв'язку граничної задачі (1) у плоскому випадку. Розв'язання задачі (2) здійс- 
нено методом розкладення за неортогональними функціями, і досліджено чисельну збіжність ітерацій. 

Досліджено похибку заміни рівняння сили тяжіння, якому реально задовольняють граничні дані зада- 
чі (1) на сфері С з радіусом Землі, на рівняння Лапласа: вона є різницею розв'язків відповідних гранич- 
них задач типу (1) для рівняння сили тяжіння та рівняння Лапласа. При розмірі площі гравіметричних 
зйомок 100 км ця похибка, як показано в |2|, варіює в межах 75 мГал 2 є 5 81 мГал за умови перерахунку 
значень сили тяжіння у вільному просторі за формулою Пуассона |31. 

Виходячи з властивостей рівняння сили тяжіння і умови Ле 20 увігнутості поверхні спостережень, 
істинне значення сили тяжіння завжди буде меншим, аніж одержане за припущення його гармонічності. 
Порівняння відповідних похибок для сфери С виявляє, що аналітична модель для сили тяжіння значно 
краще апроксимується диференціальним рівнянням типу Клейна-Гордона, ніж рівнянням Лапласа. 

Але у гравіметричній практиці при перерахунку в зовнішньому просторі розв'язують не внутрішню 
граничну задачу, а зовнішню. Максимальна величина похибки перерахунку у зовнішньому просторі щодо 
області С дорівнює тах є ,, --(а)/4, де 2(а) - значення сили тяжіння на поверхні ідеально сферич- 
ної Землі. Так, припущення гармонічності сили тяжіння в перерахунку на висоту 6000 км генерує похиб- 
ку 245 Гал. 


Завдяки розходженню векторів складових традиційних аномалій сили тяжіння дя щ Я -|8 , Їх величи- 





ни можуть значно відрізнятися від точних значень аномалій сили тяжіння. Якщо навіть аномалія сили тя- 
жіння |У -|8| високоточна, її не можна перерахувати в зовнішньому просторі: диференціальне рівняння 


сили тяжіння істотно залежить від швидкості зміни кутів між напрямом аномалії сили тяжіння і коорди- 
натними осями, а в локальних аномаліях цей напрям може значно варіювати, і всяке наближене припуще- 
ння про його поведінку істотно спотворює перерахунок аномалій сили тяжіння. Тому практично важливо, 
що вираз |--|8| співпадає з великою точністю з проекцією вектора аномалії сили тяжіння на напрям век- 


тора У. 

Чисельне моделювання. Проведено чисельне розв'язання нелінійної граничної задачі (1) для плос- 
кого і тривимірного випадку з метою перевірки збіжності ітераційного процесу (2) на модельних задачах 
з відомими точними розв'язками в горизонтальній площині 5. В процесі ітерацій застосовано метод фун- 
даментальних розв'язків в точках поверхні 5, на рівномірній сітці з кроком 10 км. 


За перше наближення напрямних косинусів Ш, узяті грубі наближення 0, 0, 1,а в іншому варіанті ро- 


зв'язку - 0.8, 0, 0.6. Кожна з граничних задач розв'язувалася для різних перших наближень і різного по- 
ложення на площині допоміжних точок, що визначають фундаментальні розв'язки. Число ітерацій у всіх 
випадках не перевищувало 10. Порівняння точних і наближених значень МГПСТ виявляє, що їх граничні 
значення у всіх точках співпадають з точністю до 107", змінюються в межах 0,07-0,133, а вони самі спів- 
падають з меншою точністю, що пояснюється поганою обумовленістю взагалі граничних задач з гранич- 
ними умовами, що містять похідні (для них на відміну від задачі Діріхле несправедливий принцип макси- 
муму). Цю обставину, не пов'язану із збіжністю ітерацій, слід враховувати при вирішенні практичних за- 
дач. 
Обчислено розв'язок внутрішньої граничної задачі 


Ли(х,) 20, х, є С, |вгайці, - 5) (3) 


для одиничного куба С, в якого зовнішня нормаль має 6 різних напрямів, гранична функція має вигляд 
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бувай початок координат (0,0,0) розташовано на відстані 0,2 від нижньої кромки куба С,а 
хуорху 5 

24 допоміжні точки узяті на кубі зі стороною |,4. За перше наближення напрямних косинусів (/, в ітера- 

ційному процесі (2) узяті сталі 0, 0 1 1. Після 5-Ї ітерації граничні значення задовольнялися достатньо 


(похибка менша 196), але функція і її похідні не мали нічого спільного з точним розв'язком 
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У с----ці---- і її похідними. Збільшення числа ітерацій до 24 не збільшило точність наближення. 
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Додали 25-й фундаментальний розв'язок, який співпадає з точним розв'язком граничної задачі (3) і після 
4-ї ітерації і граничні умови, і значення градієнтів і їх похідних задовольнялися з точністю 107. 
Збіжність ітерацій (2) окремо вивчено шляхом розв'язання зовнішньої граничної задачі (1) для оди- 


ничного круга з центром в початку координат (точний розв'язок задачі - и Іпу/х - х; , пОЧаткові набли- 
ження коефіцієнтів 7 з со(г-нФ), а з8іп(7--Ф), де г - центральний кут, Ф - міра збурення початкового 
напряму градієнта шуканої функції), з обривом ітерацій по досягненні умови я - |ргайм"| -1 «107. 


Виявлено, що для збіжності ітерацій (2) необхідно не лише, щоб точний розв'язок задовольняв умові |4| 


іп (М,І)» 0, (Ф 


хє5 
але 1 вибір такого початкового наближення, для якого в ітераціях (2) не виникали б коефіцієнти а; , ЯКІ 


не забезпечують виконання цієї умови. Умова (4) означає, що напрям зовнішньої нормалі М не має бути 
ортогональним напряму / косої похідної, тобто / в жодній точці границі 5 не має співпадати з дотичним 
напрямом. Ця умова виконується для будь-якої скінченої частини площини 5, якщо поле створене точко- 
вим джерелом, а для всієї нескінченної площини вона не виконується. 

При розв'язанні практичних задач шукають не потенціал И/ сили тяжіння, а потенціал Ї/ аномалії 
сили тяжіння, пов'язаний з ИЙ через МИ -/ --У, де У - потенціал нормального значення сили тяжіння, 
вважають відомим. За лінійних наближень граничних задач граві- і магнітометрії цей перехід до аномалії 
не змінює ні рівняння, ні вигляд оператора граничних умов, необхідно лише за граничну функцію у (5) 
узяти аномалію сили тяжіння. У нелінійних граничних задачах перехід до аномалії змінює ліву частину 
граничних умов (1) через те, що операція обчислення модуля градієнта не дистрибутивна: 
(гад МИ  егаді 2 --|вгад М| ємаа М |. 

Якщо гранична функція у (5) - модуль градієнта потенціалу Й/ сили тяжіння, для визначення поте- 
нціалу аномалії // одержуємо таку нелінійну граничну задачу 

ли(х,)- 0, х є С, ці ам ОО? НО? о У У, У 0, У, И, ) у (5)-|зтадУЇ, (5) 
граничні значення якої за нульового потенціалу нормальної сили тяжіння співпадають з граничними умо - 


вами (1). З'ясування чисельної швидкості збіжності задачі (5) представляє окремий практичний інтерес. 
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